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Precise measurement of the Higgs boson properties are important issues for the In- 
ternational Linear Collider (ILC) project to understand the particles mass generation 
mechanism which strongly related to the coupling with the Higgs boson. Large Hadron 
Collider (LHC) [T] experiments exclude the large area of the predicted Higgs mass 
region and their results indicate that Higgs boson mass will be light. Even if LHC 
discovers the Higgs like particle by the end of 2012, Higgs will be identified by the 
high precision measurement of the Higgs boson properties in ILC and also Higgs mea- 
surement verifies the correctness of standard model (SM) or gives some hints toward 
its beyond. In this study, we evaluate the measurement accuracies of Higgs branching 
fraction to the H — >• hh, cc and gg at the center-of-mass energy of 250 and 350 GeV. 

1 Introduction 

Higgs branching ratio (BR) measurement is 
one of the important issue of the Interna- 
tional Linear Collider (ILC) project, which 

is strongly related to the coupling strength Table 1: Expected SM Higgs BR at the 
with particles and reveal their mass genera- Higgs mass of 120 GeV in PYTHIA [2]. 
tion mechanism. Even the Large Hadron Col- 
lider (LHC) [1] will discover the Higgs boson 
in a few year, ILC can confirm whether that 
is the standard model (SM) predicted one or 
not, and find some hints toward its beyond. 
LHC experiment accumulate the total inte- 
grated luminosity up to 5 fb~^ by the end of 
2011, LHC gradually exclude the heavy Higgs 
mass region and indicates the light Higgs from 
the combined results of ATLAS [3] and CMS [4] (115 < Mr < 135 GeV). In this region, 
we obtain the maximum production cross section around the center-of-mass energy (^/s) of 
250 GeV and Higgs boson BR significantly varies depending on the Higgs mass, especially 
main decay channel shifts from H —> bb to WW* around the Mh = 140 GeV, as shown 
in Fig. [1] (a) and (b). In this study, we evaluate the measurement accuracies of the Higgs 
BRs of if — > bb, cc and gg with assuming the Higgs mass of 120 GeV and integrated lu- 
minosity (£) of 250 fb~^, using the International Large Detector (ILD) [S] full simulation 
at y/s ~ 250 GeV. In addition, we also consider the operation at the CM energy of 350 
GeV, which increase the contribution of W/Z fusion process but Z/H will be boosted, for 
taking into account the energy staging option in ILC project. Considering the LHC results, 
we also evaluate the Higgs BR accuracy at the indicated mass region by extrapolating the 
Mh = 120 GeV result. 
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Figure 1: (a) Higgs production cross section as a function of CM energy and (b) SM predicted 
Higgs BR as a function of Higgs mass. 



2 Signal and Backgrounds 





Figure 2: Higgs production processes categorized with Z decay channels: (a) neutrino (yvYT)^ 
(b) hadronic (qqH) and (c) leptonic (i'^i^H). 



Figure [2] shows the Higgs production diagrams and ZH analysis procedures are cate- 
gorized with the Z decay channels; ZH -> i/DH (Neutrino), qqH (Hadronic) and £^i~H 
(Leptonic) which mainly form di-jet, four-jet and di-lepton+di-jet final states, respectively. 
We also consider the W/Z-iusion process for the neutrino and leptonic channels which 
has larger contribution at -^/s = 350 GeV. Assuming the Mh = 120 GeV, we obtain 
the maximum production cross section around the ^/s = 250 GeV through the Higgs- 
strahlung (ZH) process and mainly decays to b quarks (H — > 66). In order to maximize 
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the Higgs production cross section, we employ the left-handed electron beam polarization; 
(e+, e") = (+30%, -80%), at the CM energies of both 250 and 350 GeV with the integrated 
luminosity of 250 fb^^. As background, we consider the following 2/ and 4/ final state SM 
Backgrounds: e+e^ — > W^W^ , ZZ and qq. In addition, we also take into account for the 
e+e" — ?> tt background only for the \/s = 350 GeV. 

3 Analysis Framework 

Since ZH final state forms multi-jet, thus jet clustering, jet energy resolution and quark 
flavor-tagging are crucial for the Higgs hadronic decay channels [H — >■ 66, cc and gg) analy- 
sis. In order to achieve the best jet energy resolution, ILD adopt the Particle Flow Algorithm 
(PFA) [6], which can achieve the best jet energy resolution; charged tracks energy is mea- 
sured by tracker instead of the calorimeter and only neutral particles energy is measured 
by calorimeter with avoiding the cluster overlapping and double counting. ILD detector 
design is well-suitable for the best PFA performance with adopting the finely segmented 
calorimeter, large tracker radius and strong magnetic field. For the simulation study, we use 
the ilcsoft v01_06 [7] ILC common software package. At fist we generate the MC event 
samples with the Whizard [S]. Then we perform the ILD full detector simulation with the 
Mokka [3] package assuming the ILD detector model (ILD_00). Generated hits are digitized 
and reconstructed with Marlin [1^ package and perform the PFA (PandraPFA [6]). We also 
employ the LCFIVTX [TI] flavor tagging package embedded in Marlin to identify the quark 
flavor of iJ — )• 66 and cc. For the mass production of the full simulation samples, we use the 
GRID [12] resources for the ILD LOI study [5] and flnally saved as ILC common flic format 
(LCIO [13]). 

4 Event Reconstruction and Background Reduction 

4.1 Neutrino Channel {vvH) 

4-. 1.1 Jet Reconstruction 

For the vvH channel analysis, we apply the jet reconstruction forcibly merged into the 
two-jet which comes from H — > 66, cc and gg. 

4-. 1.2 Background Reduction 

In this channel, ZZ — > i/iyqq and WW — > i^iqq will be the main backgrounds. To suppress 
these backgrounds, at flrst. we apply the missing mass cut; 80 < Mmiss < 140 GeV to 
suppress the ZZ to leptonic or hadronic decay backgrounds, since Mmiss should be consistent 
with the Z mass in Z ^ channel. Then we use the following kinematic variables cut: 
transverse momentum; 20 < Pt < 70 GeV, longitudinal momentum; \Pi\ < 60 GeV and 
maximum momentum; Pmax < 30 GeV to suppress the qq background. Number of charged 
tracks cut; Nchd > 10 well reduce the backgrounds including energetic leptons. We also 
apply the jet clustering y-value cuts; 0.2 < Y12 < 0.8 and 5^23 < 0.02, which are the 
y-value thresholds from two- to one-jet or two- to three-jets and different of number of 
jets backgrounds arc reduced with these cut. Finally we apply the di-jet mass cut; 80 < 
Mjj < 130 GeV which correspond to Higgs mass. After applying all the cuts, we apply 
the likelihood variable cut using following input variables; Mmiss 1 number of reconstructed 
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particles {Nppo), Pmax, Pi and Mjj, and we select the LR > 0.375 where we obtam the 
maximum signal significance. In order to optimize the cut positions for the analysis at 
^/s = 350 GeV, we change the following variables cut positions: 50 < Mmiss < 240 GeV, 
10 < Pt < 140 GeV, \Pi\ < 130 GeV and LR > 0.15. 



Table 2: Summary of the background reduction in viyH channel at the y/s = 250 and 
350 GeV. 



(GcV) 


250 


350 


Cut 


Sig. 


Bkg. 


Sig. 


Bkg. 


Gen. 


19360 


44827100 


26307 


20855900 


All cuts 


6731 


19058 


12338 


71918 


LR cut 


4753 


3593 


9302 


10029 


Significance (EfT.) 


52.0 (24.5%) 


66.9 (35.4%) 



4.2 Hadronic Channel (qqH) 

4-2.1 Jet Reconstruction and Pairing 

For qqH channel, we apply the four-jet reconstruction forcibly. Then we calculate following 
value to determine the Z or H jet pair candidates: 

^i^) , (1) 

where M12 and M34 are reconstructed di-jet invariant masses and ctz/h are the width of Z 
and H mass distribution. We select the minimum j^t pairs as the best candidate of Z 
and H. 




4-2.2 Background Reduction 

After the jet clustering and pairing, we apply the background reduction. At first we apply 
the cut to reduce the wrong combination pairs; x^ < 10- Then wc require the following 
cuts to suppress the leptonic events: the number of charged tracks cut; Nchd > 4 and y- 
value cut Y34 cut, which is a y- value threshold from three- to four-jet; — log 134 < 2.7. As 
event shape cuts, we employ the following variables cut: thrust and its cosine of thrust 
angle; thrust < 0.9, \ cos 9thrust\ < 0.9 and the angle between the Higgs candidate jets; 
105 < Oh < 160° to suppress the qq backgrounds. Finally we apply the di-jet mass cut M12 
and M34 and likelihood variable cut which is calculated with following input variables; thrust, 
Mz, Mh, Oh- We select the likelihood variable cut as LR > 0.2. For the ^/s = 350 GeV, we 
optimize the cut position as follows; thrust < 0.85, 70 < Oh < 120°, 80 < Mz < 100 GeV, 
105 < Mh < 130 GcV and LR > 0.1. 

4.3 Leptonic Channel i7) 

For the £^i^ H channel analysis, at first we identify the di-lepton, then we apply the di-jet 
reconstruction forcibly for remaining particles. 
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Table 3: Summary of the background reduction in qqH channel at the \/s = 250 and 
350 GeV. 



(GeV) 


250 


350 


Cut 


Sig. 


Bkg. 


Sig. 


Bkg. 


Gen. 


52507 


44827100 


36099 


21222700 


All cuts 


16350 


411785 


9447 


44400 


LR cut 


13726 


166807 


8686 


25393 


Significance (Eff.) 


32.3 (26.1%) 


47.1 (24.1%) 



4-. 3.1 Di-lepton Identification 

We apply the following di-lepton identification for electrons and muons from the differ- 
ent aspects in the energy deposition in the calorimeter; EecalI ETotai > 0.9 and 0.7 < 
Erotai/P < 1-2 for electrons and EECAL/Exotai < 0.5 and Erotai/P < 0.4 for muons, 
where Eecal, Erotai and P denote the ECAL energy associated with a track, total energy 
deposit in whole calorimeter and track momentum. 



4-. 3. 2 Background Reduction 

After the di-lepton identification, we apply the background reduction with following cut 
variables, which is summarized on the Tabled At first we apply the di-lepton mass (Mu) 
cut which should be consistent with the Z mass: 70 < M« < 110 GeV for electron and 
70 < Mu < 100 GeV for muons, respectively. Then we apply the Z fiight direction cut: 
I cos 0^1 < 0.8 to suppress the forward region backgrounds. Finally we require the di-jet 
mass {Mjj) and recoil mass (Mrec) cuts to select the Higgs candidate signal: 100 < Mjj < 
140 GeV, 70 < Mrec < 140 GeV for electron; 115 < Mjj < 140 GcV and 70 < M^ec < 
140 GeV for muon, respectively. 



Table 4: Summary of the background reduction in H channel at the -y/s = 250 and 
350 GeV. 



(GcV) 


250 


350 


Cut 


Sig. 


Bkg. 


Sig. 


Bkg. 


Gen. (e) 




3137 


4512520 


2740 


3822410 


Gen. (//) 




2917 


4512520 


1789 


3822410 


All cuts (e) 




1184 


1607 


567 


590 


All cuts (fi) 




1365 


983 


638 


465 


Significance (Eff.) 
Significance (Eff.) 


(e) 
(m) 


22.4 (37.8%) 
28.2 (46.8%) 


16.7 (20.7%) 
19.2 (35.7%) 
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5 Measurement Accuracy of Branching Fraction 



5.1 Template Fitting 

In order to evaluate the measurement accuracies of the Higgs branching fraction (BR) , we 
employ the template fitting method [M]. At first, we prepare the flavor-likeness template 
samples which is calculated from the LCFIVTX output xi.2'- 

X — likeness = (2) 

(1 - Xi)(l -X2) 

where xi^2 represents the flavor tagging output from LCFIVTX for di-jet. We assume the 
Poisson statistics {Pmk) for each bin (i, j, k) of the template samples: 

p.,k = (/i ^ E X ^ NfX") , (3) 

where Nf°j^'^ is the number of entries in (i, k) bin. represents the sum of the 

number of entries at bin (i, k) in each template sample: 

iV^^''^*^^ E r,.A^„ (4) 

s—hh,cc,gg,bkg 

where Nf-j. represents the number of entries at the (i, j, k) bin in iJ — hh, cc, gg and 
background template sample (N^^j^), which includes the number of entries of SM background 
and Higgs to none hadronic decays. Vs represents the fitting parameters of r^b, Tcc, fgg and 
fbkg , where they are the ratios of number of entries in Higgs hadronic decays of _ff — > 66, cc 
and gg after the background reduction to the entries predicted from the SM Higgs BR. r^kg 
is a normalization factor for the SM background and other Higgs none hadronic decays, 
which is fixed to be 1 from the assumption that the SM backgrounds are well understood. 
Finally we apply the template fitting with minimizing the following log-likelihood variable 
L calculated from the product of the probability Pijk in each bin: 

L = - log n P^i'' = " E ^'^■'^-) • (5) 

\i,j.k J i,j,k 

5.2 Measurement accuracies of BR 

To evaluate the measurement accuracies of the sigma times BRs for 66, cc and gg, we apply 
the 1000 times template fitting Toy-MC and obtain the fitted results of rg (s = bb, cc, gg): 

azH ■ BR{H -^s) = r, x azH^''' ■ BR{H ^ sf^' , (6) 

where azH is a Higgs production cross section, (Jzh^^^ and BR{H — > s)^^^ are cross section 
and BR of F ^ s predicted in SM. From the Eq. El we obtain the measurement accuracies 
of Higgs BR from following equation: 

Ac7BR{H s) 

aBR V V '"s y V '^ZH 

here we assume the 2.5% of cross section measurement uncertainty {lS.a z h I <^ z h) estimated 
from the recoil mass study |15j . Summary tables of the measurement accuracies of Higgs 
BR are shown in Table [5] and |6] for ^/s = 250 and 350 GeV, respectively. 
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Table 5: Summary of template fitting results and BR measurement accuracies at the -^/s = 
250 GeV. 





vDH 


qqH 


e+e-H 


l^+l^rH 


comb. 


nb 


1.00±0.016 


1.00±0.015 


1.00 ± 0.039 


1.00 ± 0.33 


1.00±0.012 


Tec 


f.00±0.12 


1.00 ±0.12 


0.98 ± 0.29 


1.01 ± 0.24 


1.00±0.09 




0.99±0.14 


1.00±0.13 


0.99 ± 0.35 


1.00 ± 0.21 


l.OOiO.lO 


aBR{bb)/a'='" (%) 


65.7±1.1 


65.7±1.0 


65.7 ± 2.6 


65.7 ± 2.2 


65.7±0.7 


(tBR{cc)/<t^" (%) 


3.59±0.43 


3.61±0.44 


3.53 ± 1.03 


3.63 ± 0.85 


3.60±0.31 


aBR{gg)/a^''' (%) 


5.46±0.76 


5.48±0.76 


5.45 ± 1.94 


5.49 ± 1.14 


5.47±0.54 


AaBR{bb)/<TBR (%) 


3.0 


2.9 


4.7 


3.3 


2.7 


AaBR{cc)/aBR (%) 


12.2 


12.3 


29.3 


23.5 


8.9 


AaBR{gg)/aBR (%) 


14.2 


14.1 


35.6 


20.7 


10.2 



Table 6: Summary of template fitting results and BR measurement accuracies at the y/s = 
350 GeV. 





wH 


qqH 


e+erH 


11+ ^rH 


comb. 


nb 


1.00±0.012 


1.00±0.015 


1.00 ± 0.056 


1.00 ± 0.051 


l.OOiO.OlO 


^cc 


l.OOiO.lO 


0.99±0.11 


1.02 ± 0.26 


1.02 ± 0.32 


1.00±0.07 




l.OOiO.lO 


1.00±0.13 


0.97 ± 0.35 


0.97 ± 0.35 


1.00±0.08 


aBR{bb)/a'='" (%) 


65.7±0.8 


65.7±1.0 


65.6 ± 3.68 


65.6 ± 3.32 


65. 7± 0.6 


aBR{cc)/a^'" {%) 


3.60±0.35 


3.68±0.26 


3.68 ± 0.94 


3.66 ± 1.16 


3.59±0.26 


cjBR{gg)/cj^^' (%) 


5.48±0.53 


5.49± 0.74 


5.32 ± 1.91 


5.35 ± 1.94 


5.48± 0.43 


AaBR{bb)/aBR (%) 


2.8 


2.9 


6.1 


5.6 


2.7 


AaBR{cc)/(jBR (%) 


10.1 


11.2 


25.6 


31.7 


7.7 


AaBR{gg)/aBR (%) 


9.9 


13.7 


36.0 


36.3 


8.2 
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6 Extrapolate to the different Higgs masses 



Since latest LHC results indicate the light Higgs (115 < Mh < 140 GeV), we should extend 
the results at the Mh = 120 GeV to the prospective mass region. In order to evaluate the 
Higgs BR measurement accuracy at other Higgs masses, we extrapolate the results at the 
Mh = 120 GeV with Eq. [7] while we assume the same background reduction efficiency for 
other mass. 

[ AaBR \ ^ f^GBR \ ^ / ai2o ■ BR{s)^ 

\ aBR ^''J \ aBR J \/ a.j^ ■ BR{s)m^ ' ^ ' 

where (Jmu and BR{s)]\j^ denote the cross section and BR of iJ — > s at Mh, as shown in 
Fig. [21(a). In addition, we also compile the H WW* Aj result [TB], even though this 
study assume the electron right-handed polarization: (e"*", e~) = (—30%, +80%) to suppress 
the e+e~ — J- W~^W~ backgrounds. Figure |3](b) shows the expected measurement accuracies 
of the Higgs BR at each mass and summarized on the Table [71 




Figure 3: Higgs mass dependence of (a) cross section times Higgs BR and (b) measurement 
accuracies of the Higgs BR extrapolated from the result at Mh = 120 GeV. 



7 Conclusion 

We evaluate the measurement accuracies of the Higgs BR for H — bb, cc and gg channels. 
With the template fitting analysis, we obtain the measurement accuracies of ct x BR for 
H — > bb, cc and gg as 3%, 9% and 10%, respectively. We also estimate the Higgs BR 
measurement accuracies at the prospective mass region by LHC with extrapolating the 
results at the Higgs mass of 120 GeV. 
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Table 7: Summary of the extrapolated accuracies of the Higgs BR from the result at the 
Mh = 120 GeV, y/s — 250 GeV to prospective masses calculated by HPROD jjj and 
HDECAY [TB]. Here 2.5% of the azH uncertainty is also included. 



Mh (GeV) 


120 


130 


140 


a(fb-^) 


354.3 


279.9 


203.1 


Modes 


BR{%) 


aBR 


AaBR 

rrRR 


BR{%) 


aBR 


AaBR 
aBR 


BR{%) 


aBR 


AaBR 

aBR 


H -^bb 


66.5 


235.6 


2.7% 


51.2 


143.3 


3.5% 


33.0 


67.1 


5.1% 


H ^ cc 


2.9 


10.4 


8.1% 


2.3 


6.3 


10.4% 


1.5 


3.0 


15.2% 


H^gg 


8.2 


29.2 


9.0% 


7.5 


21.0 


10.6% 


5.7 


11.5 


14.3% 


H W^W- 


13.6 


48.3 


15.7% 


29.4 


82.4 


10.3% 


49.2 


99.8 


9.3% 
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